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Abstract: 

1 . A high immediate mortality rate of released animals is an important cause of translocation 
failure ("release cost"). Post-release dispersal (i.e. the movements from the release site to 

25 the first breeding site) has recently been identified as another source of local translocation 
failure. In spite of their potential effects on conservation program outcomes, little is known 
about the quantitative effects of these two sources of translocation failure and their 
interactions with environmental factors and management designs. 

2. Based on long-term monitoring data of captive-bred North African houbara bustards 
30 (hereafter, houbaras) over large spatial scales, we investigated the relative effects of release 

(e.g., release group size, period of release), individual (e.g., sex and body condition) and 
meteorological (e.g., temperature and rainfall) conditions on post-release survival (N = 957 
houbaras) and dispersal (N = 436 houbaras). 

3. We found that (i) rainfall and ambient air temperature had respectively a negative and a 
35 positive effect on houbara post-release dispersal distance, (ii) in interaction with the release 

period, harsh meteorological conditions had negative impact on the survival of houbaras, 
(iii) density dependent processes influenced the pattern of departure from the release site 
and (iv) post-release dispersal distance was male-biased, as natal dispersal of wild birds 
(although the dispersal patterns and movements may be influenced by different processes 
40 in captive -bred and in wild birds). 

4. Synthesis and applications. Overall, our results demonstrate that post-release dispersal and 
mortality costs in translocated species may be mediated by meteorological factors, which in 
turn can be buffered by the release method. As the consequences of translocation programs 
on population dynamics depend primarily upon release costs and colonisation process, we 

45 suggest that their potential interactions with meteorological conditions be carefully 
addressed in future programs. 
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50 Introduction 

Although high immediate mortality of released animals is an important cause of translocation 
failure ("release cost", Tavecchia et al. 2009), post-release dispersal (an analogue of natal 
dispersal, defined here as the movements from the release site to the first breeding site) was 
identified as another source of local failure (Le Gouar et al. 2008). At the local (population) 

55 scale, translocation may serve to enhance population growth or improve population viability 
(Seddon et al. 2012). At larger scales, such as the metapopulation scale, translocation may 
promote the connectivity of local populations and/or extend the species' range (Seddon et al. 
2012). Although dispersal is a key component of metapopulation dynamics, it is also 
associated with local failure (Mihoub et al. 2011). Furthermore, if dispersal entails survival 

60 costs (Bonte et al. 2012), it may partially explain patterns of mortality associated with 
restoration projects (Tweed et al. 2003). Overall, the factors affecting mortality and dispersal 
following release events strongly influence the success or failure of population management 
and conservation strategies (Le Gouar et al. 2008, Mihoub et al. 201 1, Oro et al. 201 1). 
While adequate post-release monitoring are now prevalent (Seddon et al. 2007), temporal and 

65 spatial replicates of release events are still uncommon (Le Gouar et al. 2008). Consequently, 
conservationists rely primarily on short-term survival and dispersal data to assess the effects of 
release on population dynamics and to inform release strategies. However, the short-term 
survival of a released population might poorly reflect persistence over the long-term 
(Armstrong et al. 1999). Moreover, individual decisions and abilities along the dispersal 

70 process are differently affected by proximate factors (Clobert et al. 2009). In particular, if 
variations in dispersal patterns and movements are shaped by the individual phenotype and 
environmental conditions, this may result in different dispersal patterns and movements in 
captive-bred and wild individuals, having in turn an impact on translocation success (Clobert et 
al. 2009). Hereafter, we will refer to the three stages of dispersal following Bonte et al. 



4 



Downloaded from http://biorxiv.org/on September 18, 2014 



75 (2012)'s general definitions and nomenclature for wild populations (i.e., departure from natal 
area, transfer -the movement per se and settlement at the reproduction site), although adapted 
to the case of translocations (i.e., considering departure from the release site, see Le Gouar et 
al. 2012). Assessment of the multiple factors affecting survival and dispersal over different 
temporal and spatial scales is necessary for conservationists to comply with the objectives of 

80 the program (e.g. population persistence or connectivity among populations) and to determine 
the sustainability of translocation over time (Armstrong & Seddon 2008). 
Post-release survival and dispersal of individuals may be influenced by two main types of 
factors: 1) those relating to the initial conditions of the release (for instance and non- 
exhaustively: release group size, period of release; hereafter, release factors) and 2) those 

85 relating to individual characteristics (hereafter, individual factors), such as sex, age, and body 
condition. Additionally, although it is likely that local meteorological factors can also greatly 
influence the survival and dispersal of translocated animals, they have, to our knowledge, 
never been investigated. 

This paper evaluates the importance of short- and long-term processes and the relative effects 
90 of release, individual and local meteorological factors on the survival and dispersal of a 
captive-bred species in a reinforcement (i.e. an addition of individuals to an existing 
population, IUCN 2012) conservation program. We benefited from an important conservation 
and release effort involving long-term and large spatial-scale monitoring of captive-bred North 
African houbara bustards (Chlamydotis undulata undulata, Family Otididae; hereafter, 
95 houbaras) in Morocco. 

In numerous species, (i) meteorological factors strongly influence individual survival (Nevoux 
et al. 2008), particularly in low quality, weak or inexperienced individuals (Robert et al. 2012); 
and (ii) released individuals suffer high immediate mortality after release (Tavecchia et al. 
2009). Based on these two empirical findings, we predict that survival of released houbaras 
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100 will be negatively affected by harsh meteorological factors (such as abundant rainfall and/or 
low temperature) and that their effects will be strongest shortly after release due to the 
inexperience of houbaras in their release habitat. 

Although most theoretical predictions deal with dispersal rate but not with the timing of 
dispersal or the dispersal distance, most classical hypotheses on dispersal rates in patchy 

105 habitats were generalized to movement patterns and can be applied to non-colonial species and 
continuous habitat (e.g. Debarre & Gandon 2010). Accordingly, we predict that the time 
elapsed between release and departure will be negatively correlated with release group size, a 
surrogate of the local density at the time of release, as (i) the houbara's release protocol 
involves very large groups and (ii) intraspecific competition has been identified as an 

110 important ultimate cause of dispersal in both theoretical works and empirical studies of many 
species (Perrin & Mazalov 2000), including houbara bustards (Hardouin et al. 2012). In 
addition, recent evidence suggests that natal dispersal distance is male-biased in this species 
(Hardouin et al. 2012). Therefore, if sex-specific dispersal behaviour is not altered by captive- 
breeding or release protocols and conditions, we predict that post-release dispersal distance of 

115 houbaras will be male-biased. 

Methods 

GENERAL METHODS 
Study sites and biological model 
120 The North African houbara bustard is a non-migratory, sexually dimorphic species, with males 
being larger and heavier than females. It exhibits an exploded-lek mating system, where 
display males are less tightly aggregated than in a classical lek (Hingrat et al. 2008). This 
vulnerable species (BirdLife International 2012) has suffered severe declines in wild 
populations over the last 40 years. Conservation concern led in 1995 to the establishment of 
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125 the Emirates Centre for Wildlife Propagation (ECWP) in Morocco (Lacroix 2003). A 
comprehensive strategy combining ecological research, reinforcement using captive-bred 
houbaras and hunting management was established in order to restore houbara populations and 
allow regulated traditional Arab falconry. A partnership was established between ECWP and 
Moroccan authorities for the conservation and management of houbara populations over an 

130 area of 74 400 km 2 (Fig. 1). This study area represents between 8 and 13% of the species 
distribution range in North Africa (using Species Distribution Modelling with houbara 
presence thresholds of 0.4 and 0.5, A.- C. Monnet, unpublished data). Hunting was banned 
between 2000 and spring 2005 and then restricted to the wintering period (i.e. October to 
January) in delimited hunted areas totalizing 60% of the managed area. The present study was 

135 conducted from March 2001 to February 2010, then comprising a non-hunted period (2001 to 
spring 2005) and a hunted period (autumn 2005 to 2010). 
Breeding stations, release and radio-tracking procedure 

Two ECWP breeding stations were established: one in Missour in 1995 and one in Enjil (Fig. 

1) in 2005. Between 1998 and 2012, the ECWP has released 61,555 houbaras in North Africa, 
140 primarily in Morocco (see Lesobre et al. 2008 for insemination and rearing protocols, housing 

and pre-release conditions). Houbaras were released during two periods: from September to 

November (hereafter, "autumn release") and from March to May (hereafter, "spring release"). 

In autumn, houbaras are mainly released in non-hunted areas to avoid direct mortality and 

disturbance due to hunting. The releases varied by sites (Fig. 1) and group size (from 10 + 3 to 
145 292 ±101 houbaras per release) over time. Before release, each individual was weighed (± 1 

g), blood-sampled (for molecular sexing) and tagged with a metal ring or subcutaneous 

electronic micro-chip (TROVAN LID 100, DorsetID). 

One thousand and four hundred houbaras (out of the 20,652 houbaras released during the study 
period) were fitted with necklace battery-powered transmitters with mortality signals (llg 
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150 model RI-2B-M, 20g model RI-2D-M, Holohil System Ltd.). Sixty-five birds were recaptured 
before their batteries died and fitted with new 20g transmitters. Houbara bustards were 
monitored from the ground and by aerial telemetry (details in Hardouin et al. 2012). Birds were 
primarily ground-based monitored at their last location within a radius buffer of 5km. When 
not detected, birds were searched for using aerial telemetry (detection up to 80km, 303 ±138 

155 flying hours per year) for six months following an approximate radius of 40km from its last 
location and were systematically searched for during other aerial telemetry procedures (i.e. 
when searching for other "lost" birds). After six months without any location, a bird is 
considered as "lost" and searching is ceased. Locations were taken with GPS (accuracy + 5 
meters for terrestrial telemetry and ± 200 meters for aerial telemetry). 

160 Two standard monitoring protocols were applied: 1) from their release to 90 days (±10 days) 
post-release, houbaras were searched for twice a week over the entire study area in Morocco 
(Fig. 1), and 2) after 90 days, houbaras were searched for once a week. The first protocol was 
primarily used to get accurate measures of movement to analyse the departure and transfer 
phases of the short-term post-release dispersal (see below). Both protocols were used in 

165 combination to determine the distance from the release site to the first breeding site (i.e. 
settlement distance) and to build capture histories for survival analyses. For survival analysis, 
we considered only individuals released in north-eastern Morocco (Fig. 1), as release 
monitoring in the remaining study areas ceased after 90 days. 

170 POST-RELEASE DISPERSAL 

Definition of movement parameters and transfer phase 

Short-term post-release dispersal behaviour, defined as departure and transfer within three 
months of release of houbaras, was described with two distinct metrics. 
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Metric 1 corresponds to the linear distance between an individual's release location and any of 
175 its successive locations; it measures the propensity to disperse away from the release site 
(Hardouin et al. 2012). Repeated measures of this metric over the whole monitoring period 
characterized the net dispersal distance (NDD). Metric 2 is the distance between two 
successive locations. From this metric, we calculated cumulative distance (CD), as the full path 
distances covered by an individual from the release site to the last location (Hardouin et al. 
180 2012). This metric was used to estimate departure date (see below). Distances were measured 
using ArcGIS 9.2 (Environmental Systems Research Institute) and Hawth's Analysis Tools 
(Beyer 2004). Only birds located at least 10 times and at least once per 30 days were 
considered for analysis. Birds that died or were lost during this period were excluded from 
analysis (Hardouin et al. 2012). Therefore, 436 individuals and 8,679 localizations were 
185 considered in this analysis, yielding 20 ± 8 (Mean ± SD) localizations, 92 + 8 days of 
monitoring per individual and 4.7 ± 4.3 days between localizations. 
Departure phase 

The date of departure from the release site was estimated using cumulative distances (CD). 
Departure date was defined as the date at which the distance from the release sites begins to 
190 increase with time, inducing a break on the regression line between cumulative distances and 
time (detailed procedure in Appendix SI). 

Settlement distance 

Eighty six houbaras (out of the 436 individuals) were continuously monitored from their 
release until their first breeding attempt. Their settlement distance was estimated as the linear 
195 distance from the release site to the first nest (for females) or display site (for males). Released 
females initiated reproduction at 1.6 + 0.5 years, whereas released males started displaying at 
2.1 + 0.8 years (data from our study sample), making them difficult to monitor over the long- 
term. We monitored 69 females and 17 males with respect to settlement distance. 
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Local meteorological data 

200 In the north-eastern region of the study area, a network of 10 weather stations was installed 
close to the primary release sites (Fig. SI). Ambient air temperature (°C) and rainfall (mm) 
were recorded daily. For each bird's location, we calculated the mean ambient air temperature 
and sum rainfall for the + 3 days surrounding the localization event, using data from the 
nearest weather station. Rainfall and temperature were negatively correlated (r = -0.13, P < 

205 0.001) with high temperature associated with dryness. We assessed the effects of local 
meteorological factors on houbaras during transfer using only individuals for which release site 
and subsequent localisations were near the weather stations to properly assess the effects of the 
meteorological variables (N = 223 out of 436 individuals). The mean distance between bird's 
locations (N = 3,628 localizations) and weather stations was 17.8 + 12.7 km; and did not 

210 exceed 65 km. 

Released bird body condition and age 

The 436 houbaras considered in the analysis were released at an average age of 234 + 89 days, 
weighing 1,092 + 140 g (females) and 1,482 + 235 g (males) (see Table SI for details). An 
index of body condition was extracted for each individual using the residuals of the regression 
215 of the body weight and release age (as body weight linearly increased with age at the time of 
the release; inflexion point of the growth curve = 44 + 8 days for males and 38 + 6 days for 
females) with sex as a covariate (F 2, 433 = 322.7, Adj.R 2 = 0.6, P < 0.001, Hardouin et al. 
2012). 

Statistical analysis 

220 We first investigated the effects of release and individual factors on departure date using linear 
mixed effect models. The model tested for effects of body condition (and its quadratic term), 
sex, release age, release group size and first order interactions on departure date. 'Year', 
'release site' and 'period' were added as random terms. We then tested the effects of body 
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condition (and its quadratic term), sex, release age, release group size, time after release (i.e., 

225 based on the dates of the localization events) and first order interactions on the NDDs to 
analyse the pattern of movement along the short-term post-release dispersal. We considered 
several random terms: 'individual ID' (nested into 'release site', N = 22), 'year' and 'period of 
release'. Third, we tested for effects of local meteorological factors (i.e., rainfall and ambient 
air temperature), sex, body condition and their first order interactions on NDDs on individuals 

230 for which accurate meteorological data were available (N = 223 individuals, see above) to 
analyse their effects on the short-term dispersal distance between each individual's 
localizations. The random slope terms were added as in the previous model and 'time after 
release' was considered as a random intercept to account for the repeated measures over time. 
Finally, we investigated whether settlement distance was explained by body condition at the 

235 time of the release, sex or age at first reproduction (as longer-distance dispersers might 
reproduce later). 'Year of reproduction' was added as a random slope term. For all models, we 
started from the global model (all explanatory variables and first order interactions) and 
compared its performance with submodels from which non-significant terms were deleted one 
at a time (see details in Appendix S2). NDDs and departure dates were log-transformed and the 

240 cumulative distances were square -root-transformed to meet normality and homoscedasticity 
assumptions. The freeware R 2.10.1 (R Development Core Team 2009) and the libraries lme4 
(Bates & Maechler 2009), languageR (Baayen 2009) and segmented (Muggeo 2008) were used 
for statistical analyses. 

245 SURVIVAL ANALYSIS 

We used multi-event capture-recapture modelling (Pradel 2005) using a mixture of alive 
recaptures and dead recoveries to estimate survival of radio-marked individuals (N = 957 
individuals) from March 2001 to February 2010. Capture occasions were defined as a time- 
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step of three months (N = 36 occasions). We studied the effect of sex and release period 
250 (autumn and spring) on short- (next occasion) and long-term (following occasions) survival 
probability according to years. The time dependent model was coded to possibly obtain (when 
relevant): 1) the survival in the third months after release (i.e. short-term survival) for each 
release cohort, sex and year and 2) the annual survival on the long-term (after the third months 
post-release) for each release cohort and sex. We also assessed the influence of individual and 
255 temporal covariates on survivorship. Each bird was monitored for an average of 351 + 400 
days (max = 2,721 days). Details of the sample used in this analysis and of the model structure 
are given in Table S2 & Appendix S3. 

Briefly, we considered two states to code whether an individual is alive or dead and 3 events to 
code for the observed fate of an individual at each occasion (0: not observed, 1: alive 

260 individual, 2: recovery of a dead individual). We verified the fit of the general, time-dependent 
model with program U-CARE version 2.3.2 (Choquet et al. 2009a, details in Appendix S3). A 
step-down model selection was performed using program E-SURGE version 1.7.1 (Choquet et 
al 2009b) based on the Akaike Information Criterion corrected for sample size and 
overdispersion (QAICc). When AQAICc was smaller than 2, we selected the model with the 

265 smallest number of parameters, following the parsimony principle (Lebreton et al. 1992). 
Starting from the time -dependent model (Appendix S4), we thus proceeded of three successive 
steps of model selection: 1) we performed model selection for recovery probabilities (Models 
El to E10 in Table S3); 2) recapture probabilities (Models CI to C4 in Table S4); 3) and 
survival probabilities (Models SI to S16 in Table S5). Unless indicated otherwise, estimates of 

270 survival, detection and recovery probabilities are presented as the value + standard error (S.E.). 
Selection of individual and meteorological covariates and modelling 

To test whether variation in meteorological conditions explains a significant part of the 
temporal variation in survival rate, we integrated meteorological information as temporal 
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covariates into models. At the spatial scale of our study (Figs. 1 & SI), no data on temporal 
275 variation in the abundance of houbara food resources were available. However, the availability 
of plant and animal species in the houbara' s diet is known to vary considerably over space and 
time (Bourass et al. 2012), and meteorological parameters have been used to characterize the 
environmental fluctuations faced by individuals (Nevoux et al. 2008). We used the daily 
temperature and rainfall data from the 10 weather stations (Fig. SI) to create annual 
280 meteorological indices. We then calculated the annual mean daily temperature (ADT) and 
rainfall amount (ADR) over the ten stations. After 6 years of reinforcement, regulated hunting 
on houbaras occurred from 2005 to 2009 in the winter (non-breeding) season. Although 
hunting was not our focus and autumn releases mainly occurred in non-hunted areas, we 
accounted for hunting in our analyses by integrating a temporal covariate describing the 
285 presence (1) or absence (0) of hunting in a given year in our model with meteorological 
covariates. 

We first assessed the effects of meteorological indices on the selected time-dependent model to 
estimate which one (ADT, ADR or their interaction) had the greatest effect on temporal 
variation in survival (Models Tl to Til in Table S6). The contribution of temporal covariates 
290 to the model was assessed using an analysis of deviance (ANODEV), and the amount of 
temporal variation in the focal rate explained by the covariate was evaluated by R 2 dev (Skalski 
et al. 1993). 

We also tested for the effects of individual covariates describing the initial conditions of birds 
and group size upon release. We used an index of body condition (BC) (computed as above), 
295 the age of individuals at release (Ra) and the release group size (RGr) to assess whether these 
covariates influence survival rates. We individually tested each survival parameter of the time- 
dependent selected model, resulting in 1 6 candidate models with individual covariates (Models 
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1 to 16 in Table S7). A covariate is important when the associated 95% confidence interval 
(CI) of the corresponding slope does not include zero. 

300 

RESULTS 
DISPERSAL 

Effects of release and individual factors on departure date 

Departure dates were calculated for 396 of the 436 individuals included in this analysis and 

305 occurred on average 29.7 ± 20.2 days after release (median = 24.7, range = 2-90 days). From 
the selected model (R 2 = 0.13, x 2 miii = 8.2, df =3, P = 0.04), departure date was weakly but 
significantly explained by a positive interaction between body condition and release group size 
(MCMCmean= 1.04, HPD95 {1.001 - 1.08}, P = 0.04): individuals in poor body condition 
released in a large group tended to leave the release site first. 

310 Effects of release and individual factors during transfer 

The most parsimonious model of NDDs (R 2 = 0.65, x 2 nu ii= 2387.8, df = 7, P < 0.001) revealed 
significant effects of body condition, release age and time after release and first order 
interactions of time after release with body condition, release age and sex (Table 1). 
Individuals in better body condition moved farther distances than those in poorer body 

315 condition (Fig. 2). Older released birds travelled farther distances shortly after release, but not 
later (Table 1, Fig. S2). Moreover, the pattern of movement between males and females 
differed over time (sex * time after release interaction, Table 1). Males moved shorter 
distances than females shortly after release, but over time, the pattern reversed. This interaction 
suggests an increase of NDD in males or a decrease of NDD in females. Release group size did 

320 not affect NDDs (Release group: MCMCmean= 0.0007, HPD95 { -0.0006 - 0.002 } , P = 0.4). 
Effects of meteorological factors during transfer 
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According to the selected model (R 2 = 0.76, % 2 nu ii = 16.1, df = 3, P = 0.001), body condition 
was still positively associated with movements when including meteorological factors as 
explanatory variables (estimate + SE = 0.0012 ± 0.0004, t-value = 3.1). Along monitoring, 
325 individuals covered greater distance when facing high ambient air temperature (estimate ± SE 
= 0.0133 ± 0.0045, t-value = 2.9) and low rainfall conditions (estimate ± SE = -0.0226 + 
0.0067, t-value = -3.3). 

Effects of release and individual factors on settlement distances 

The selected model (R 2 = 0.21, df = 1, x 2 nuii = 7.15, P = 0.007) revealed that settlement distance 
330 was sex-biased, with males settling farther from the release site than females (MCMCmean = 
0.61, HPD95 {0.15; 1.05}, P = 0.008). Neither body condition nor age at first reproduction 
significantly affected settlement distance. 
SURVIVAL 

Model selection procedures are detailed in Appendix S4, Tables S3 & S4 for encounter 
335 probabilities and Table S5 for survival. 
Modelling variation in survival rates 

Ultimately, the model [S aVtperiodVtt+aVtperiod2 +a2*t\ assuming short- vs. long-term- (i.e., al 
and a2, respectively), time- and release period-dependent survival was selected (Model SI in 
Table S5). In particular, short-term survival differed according to the period of release (i.e., 

340 autumn vs. spring, period 1 vs. period 2, respectively), with time-dependence on short-term 
survival for individuals released in autumn ([S al * periodl * t ] range between 0.19 ± 0.07 in 2009 
and 0.92 ± 0.06 in 2006, Fig. 3a), while short-term survival was constant over years for 
releases occurring in spring ([S al * period2 ] = 0.86 ± 0.02). Long-term survival varied according 
to year ([S a2 *tL Fig- 3b), with a relatively low survival rate in 2008 (0.72 ± 0.03). 

345 Effect of the meteorological covariates on survivorship 



15 



Downloaded from http://biorxiv.org/on September 18, 2014 



Meteorological covariates were incorporated into the selected time-dependent model 
[Sai* P erioai*t+ai*period2+a2*tl ( Table S6). The interaction between annual temperature (ADT) 
and rainfall (ADR) explained a significant part of the short-term survival for autumn releases 
([Sai* P eriodi*t*(ADT*ADR+hunting)], Fl g- 3a), and an additive effect of ADT and ADR explained 

350 a significant part of the long-term survival ([S a2 *t*(ADT+ADR+hunting)]> Fig- 3b, Models Tl in 
Table S6). The selected meteorological indices explained 83% of the deviance in the temporal 
variation of survival probabilities (F-Anodev 19 = 44.6, P < 0.001). It is important to note that 
when hunting is removed from the previous model, 76% of the deviance was still explained by 
meteorological covariates (F-Anodev 1 n = 35.7, P < 0.001). 

355 Effect of the individual covariates on survivorship 

Individual covariates were incorporated into the selected time-dependent model (Table S7). 
The models with 'Release age' [cov Ra ] as covariates (Models 1, 2 and 3 in Table S7) on long- 
term survival were preferred over all of the other models, with age at release being negatively 
related to the long-term survival rates (\S aUperiodUt+aUperiod2+a2 , t+a2<i+COVRa) ]: slope = - 

360 0.20 ± 0.07; CI {-0.05; -0.35 }, Model 1 in Table S7). 

DISCUSSION 

We found that (1) density-dependent processes (e.g. intra-specific competition) partly shaped 
the pattern of (early) post-release dispersal, as suggested by the interaction between initial 
365 density (release group size) and body condition; (2) post-release dispersal was male -biased, as 
in natal dispersal of wild houbaras (Hardouin et al. 2012); nevertheless, our results suggest that 
the dispersal decision may be influenced by different processes in wild and captive-bred birds; 
and (3) meteorological factors had critical effects on houbara survival, although these effects 
can be mitigated by appropriate timing / method of release; 
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370 Effects of the meteorological, release and individual factors on post-release dispersal 
phases 

Our results highlight the importance of meteorological conditions during the transfer phase of 
post-release dispersal and their potential impact on the colonisation processes of translocated 
populations. We found that rainfall and ambient air temperature had a negative and a positive 

375 effect, respectively, on net dispersal distance, independently of the release period. These 
meteorological factors can have indirect effects mediated through food availability. Recent 
findings showed that the houbara plant diet is dominated by young green leaves of annual 
plants, whereas ants dominate the animal diet of the species (Bourass et al. 2012). The 
availability of these food resources is driven by ambient air temperature and rainfall. For 

380 instance, while rainfall induces instant new sprout and high ant foraging behaviour, warm 
conditions are associated with low plant availability and low ant activity (e.g. Bernstein 1975, 
Whitford 1978, Lopez et al. 1993). Meteorological conditions can also directly affect flight 
ability, having in turn an impact on dispersal distances (Walls et al. 2005). 
The release group size interacted with body condition and affected departure from the release 

385 site. This finding suggests that individuals with poor body condition leave the vicinity of the 
release site sooner in larger release groups. Birds in good condition may be less affected by 
density-dependent processes (such as intraspecific competition) before their departure than 
birds in poor body condition (Matthyssen 2005). These birds may remain at the release site for 
longer to limit the release cost and may gain sufficient reserve and/or experience in their novel 

390 environment. This behaviour may increase their chance to survive, as opposed to the birds in 
poor body condition that seem to leave the release site prematurely in this high-density context. 
Previous results on wild juvenile houbaras indicate that males in poor body condition move 
farther than those in good condition during natal dispersal; distance being independent of body 
condition in wild females (Hardouin et al. 2012). Here, captive-bred houbaras show the 
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395 opposite pattern suggesting that selective pressures exerted by captivity (Robert 2009) and/or 
the release protocol (Le Gouar et al. 2012) may interact with the ultimate factors shaping 
individual dispersal strategies (e.g., inbreeding avoidance and intra-specific competition). 
Finally, in contrast to wild juveniles that already occupy favourable breeding habitats, released 
individuals may have to search more extensively for lekking or nest sites. Thus, only 

400 individuals in good condition may afford the distance-dependent cost of dispersal (Bonte et al. 
2012). Hence, older released birds had longer net dispersal distances than younger ones just 
after release, but not later. Presumably, the longer distances initially covered by older 
individuals reflect differences in muscular development or hormonal activity such as an 
increase in corticosterone with age during captivity (Dufty & Belthoff 2001). 

405 Finally, we detected male-biased dispersal in houbaras, consistent with studies of natal 
dispersal in wild houbaras (Hardouin et al. 2012). Sex-biased dispersal is generally interpreted 
in the theoretical frame of the inbreeding avoidance hypothesis (Szulkin et al. 2013). Recent 
findings in Morocco demonstrated a lack of kin association within and between leks in male 
houbara but a higher genetic relatedness in males than in females (Lesobre 2008). This finding 

410 suggests greater philopatry in males (Lesobre 2008), in contrast to our results. However, rare 
long-distance dispersal events in females (3 females dispersed more than 100 km in our study) 
suggest that they may be vectors of gene flow at the metapopulation level (Ibrahim et al. 
1996). This hypothesis is also supported by our data on breeding dispersal distance (i.e. the 
distance between two successive breeding sites): males, once settled, remain near their display 

415 sites (mean breeding dispersal distance ± SD = 2.3 ± 4.6 km; median = 438 m; N = 16 males, 
unpublished data), whereas females exhibited farther breeding dispersal distances (mean = 16.4 
± 30.4 km, median = 2.9 km; N = 39 females, unpublished data). 

Effects of the meteorological, release and individual factors on short- and long-term 
survival 
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420 Meteorological conditions strongly influence individual survival (e.g., Grosbois et al. 2008, 
Nevoux et al. 2008) and, consequently, population dynamics and viability (e.g., little bustards 
Tetrax tetrax, Morales et al. 2005). While several studies have reported substantial short-term 
mortality costs in translocated populations (Tavecchia et al. 2009), these costs have primarily 
been attributed to release methods, individual characteristics and habitat requirements. Here, 

425 we report a strong interaction between release strategy (autumn vs. spring releases) and 
meteorological conditions on short-term survival, which suggests that release costs are likely to 
be modulated by environmental stochasticity (Nicoll et al. 2003). Indeed, short-term survival 
was strongly associated with annual rainfall, ambient air temperature and their interaction, but 
only for individuals released in autumn. In particular, low temperatures in association with 

430 either flooding or aridity negatively affected short-term survival. In contrast, the short-term 
survival of individuals released in spring remained high and constant over years. Finally, we 
did not observe any difference in long-term survival between release cohorts. Annual rainfall 
and temperature also affected long-term survival, particularly in 2008, when severe flooding 
occurred in the study area. Although environmental stochasticity is in essence unpredictable, 

435 knowledge of its magnitude and its complex effects on the life histories of released individuals 
should greatly improve the accuracy of demographic projection models designed to examine 
the feasibility of translocation protocols. We thus advocate that environmental fluctuations be 
considered in future release strategies in houbaras and other translocation efforts (Robert et al. 
2007, Saether et al. 2007). 

440 Our results also showed that younger released individuals are more tolerant later in life to 
meteorological variation. The choice of the release age in translocation program is particularly 
important and conservationists are often divided on favouring either age or habitat related 
advantages in captive -bred individuals (Sarrazin and Legendre 2000). In our context, all 
individuals might pay the cost of release, as no learning relevant to survival in the wild 
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445 environment (e.g., identification of resource locations, anti-predator behaviour) can occur in 
pre-release aviaries. However, older released individuals might enjoy short-term advantages 
(e.g., in terms of competitive ability) and/or gain more experience with ageing that fade outside 
of the captive environment. Therefore, selection on older individuals might be expressed later, 
potentially explaining their higher mortality over the long-term. Time spent in captivity can 

450 alter the behaviour and/or physiology of individuals (e.g., Pedersen & Jeppesen 1990) in ways 
that negatively affect survival once released (Bremmer-Harrison et al. 2004). Further 
investigations are needed to understand the behavioural and/or physiological basis of age 
variation in survival. Finally, survival rates of captive-bred houbaras appear to approach or 
even exceed those of wild houbaras, although a survival analysis on the same time-scale with 

455 both wild and captive-bred houbaras would be necessary to fully evaluate this. Within three 
months of dispersal, wild juveniles (particularly females) exhibited overall lower survival rates 
(from 2006-2009: males = 0.86 ± 0.05 and females = 0.58 + 0.06; Hardouin et al. 2012) than 
adults, and adult survival rates (from 2006-2009: males = 0.91 ± 0.02 and females = 0.83 ± 
0.03) were similar those reported here. 

460 Conservation implications and conclusions 

In conservation translocations, post-release dispersal rate can be particularly high after release 
(Hester et al. 2008). Thus, understanding the factors influencing post-release dispersal is 
crucial in managing the trade-off between site fidelity and adaptive dispersal (Le Gouar et al. 
2012). Whereas high post-release dispersal rates can be associated with dispersal related 

465 project failures (e.g. dispersal of released individuals outside the focal populations, Le Gouar 
et al. 2008), the absence of dispersal can limit habitat choice processes or be detrimental at 
larger (e.g. metapopulation) scale (Trakhtenbrot et al. 2005). Our results suggest that short- 
term post-release dispersal can be greatly influenced by captive breeding and release factors 
(i.e. dispersal patterns are distinct from those of wild-born individuals and post-release 



20 



Downloaded from http://biorxiv.org/on September 18, 2014 



470 dispersal is influenced by timing and local release conditions) that however do not impact 
longer-term dispersal in terms of distance. Therefore, the assessment of the post-release 
dispersal at different temporal scales seems essential in translocated populations. While initial 
conditions lead to different responses between wild and captive-bred houbaras, these do not 
result in strong differences in natal and post-release dispersal distances (natal dispersal 

475 distance: wild male 19.4 + 12.8 km, wild female = 8.6 + 6.5 km; Hardouin et al. 2012). 
Moreover, sex-biased dispersal that constitutes an essential element in the connectivity and 
gene flow between populations occurs in both captive-bred and wild populations. Our results 
suggest that both the release group size and the period of release have to be carefully 
determined to limit premature departures of weaker individuals and long distance dispersal that 

480 may be induced by harsh meteorological conditions. Previous works have shown that local 
densities at the release site can also influence dispersal rate and pattern (Le Gouar et al. 2008). 
Further investigations on the effect of the local density could be helpful in refining release 
strategies. Overall, our results on the effects of the factors (including sex) influencing post- 
release dispersal distance enable general predictions of movement among individuals in 

485 translocated and natural populations. It is indeed important to account for post-release dispersal 
in population projections (Armstrong & Reynolds 2012) and to distinguish mortality from 
dispersal to adopt appropriate management strategies and maximise future establishment 
success (Tweed et al. 2003). 

Our results also highlight the importance of long-term monitoring program to assess survival 
490 cost, as short-term survival rates do not accurately reflect long-term survival. Our findings 
further suggest a strong short-term sensitivity of individuals to unpredictable meteorological 
variation in autumn, in contrast to individuals released in spring, when meteorological 
conditions are generally milder. No subsequent effect of release period was found on long-term 
survival, suggesting that once individuals gathered experience in their novel environment, 
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495 individuals from both release periods have similar survival probabilities. As the consequences 
of translocation programs on population dynamics depend primarily upon release costs, we 
suggest that their potential interactions with meteorological conditions be carefully addressed 
in future programs. Finally, and although survival and dispersal constitute critical components 
of population dynamics, recruitment and fertility rates need to be thoroughly assessed in future 

500 work to have a complete view of the demography of released houbaras. The Population 
Viability Analysis framework (Beissinger & McCullough 2002) represents the most useful 
approach to integrate the effects of meteorological fluctuations, hunting and management 
actions on the whole species life cycle. Such integration will provide reliable insights on the 
viability of the released population, and in turn, on the contribution of the conservation 

505 translocation program to the viability of the whole species. 
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Table 1: Effects of release and individual factors on net dispersal distance in houbaras. The results of the selected linear mixed-effects model are 
given, with NDD as the response variable and body condition, sex, release age, time after release and their interactions as explanatory variables. 
Parameter estimates (i.e., fixed-effects estimates), HPD intervals, MCMCmeans and p-values based on the posterior distribution (pMCMC) for the 
fixed-effects table are given. For random effects, the estimates of the residual standard deviation, the standard deviation of the random effects 
(intercept), HPD intervals and MCMCmeans are provided. 'Time' refers to 'time after release'. 



Fixed terms 


Estimate 


MCMCmean 


HRD951ower 


HPD95upper 


pMCMC 


Intercept 


1980.24 


1923.25 


391.00 


10627.22 


<.001 


Body condition 


1.17 


1.17 


1.13 


1.22 


<.001 


Sex 


0.93 


0.92 


0.82 


1.04 


0.19 


Release age 


1.18 


1.20 


1.06 


1.35 


0.002 


Time 


5.83 


5.84 


4.85 


6.94 


<.001 


Body condition x Time 


0.9994 


0.9994 


0.9991 


0.9998 


0.001 


Sex x Time 


1.1339 


1.1385 


1.0167 


1.2798 


0.008 


Release age x Time 


0.9986 


0.9986 


0.9979 


0.9992 


<.001 


Random terms 


Std.Dev. 


MCMCmean 


HRD951ower 


HRD95upper 




Individual : Release site 


0.69 


0.49 


0.46 


0.52 




Release site 


0.36 


0.33 


0.23 


0.45 




Year 


0.26 


0.29 


0.15 


0.44 




Period of release 


0.15 


0.76 


0.01 


1.57 




Residual 


0.75 


0.77 


0.76 


0.78 
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690 Figure 1: Locations of main release sites in Morocco. Houbaras were bred at two breeding 
stations: Missour (grey circle) and Enjil (white circle). 

Figure 2: Partial plot representing the interaction between time and body condition (BC, given 
as quantiles) according to the net dispersal distance (in km) of houbaras (log-transformed data 
were back-transformed). 

695 Figure 3: Annual variation in houbara survival rates: a) Survival rates estimated for the first 
time-class (i.e., the third month) and released in autumn, b) Survival rates for the second time- 
class (i.e., more than 3 months). For both, we present survival rates estimated from the time- 
dependent model and the 95% confidence intervals (white square). Survival rates (black lines) 
and 95% confident intervals (dotted lines) are estimated from the temporal covariates for a) 

700 with the interaction between ADT and ADR. [$ a i* pe riodi*t*(ADT*ADR)l an d f° r b) survival rates 
and 95% confident interval estimated with additive effect of ADT and ADR [S a 2*t*(ADT+ADR)]- 
c) Annual variation of the ambient air temperature (ADT, straight black line) and the daily 
rainfall (ADR, straight dashed line). Note that hunting stopped from winter 2001 to 2004 and 
started from winter 2005. The presence (1) / absence (0) of hunting were integrated in all the 

705 models as temporal covariate. Note that [S al » per i 0C i2] is constant = 0.86 ± 0.02 and not shown. 



31 




32 



Downloaded from http://biorxiv.org/on September 18, 2014 



o 



O — 



Body condition values 

= Lower BC 

= Lower than average BC 

= Average BC 
- ■ - ■ = Better than average BC 

= Better BC 



/ 




o - r 
0 



20 



40 



60 



~~ r 

80 



100 



Time after release (days) 



33 



Downloaded from http://biorxiv.org/on September 18, 2014 




